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ASYMPTOTIC ESTIMATE OF EIGENVALUES OF 
PSEUDO-DIFFERENTIAL OPERATORS IN AN INTERVAL 

KAMIL KALETA, MATEUSZ KWASNICKI, JACEK MALECKI 


Abstract. We prove a two-term Weyl-type asymptotic law, with error term O(^), for 
the eigenvalues of the operator ip(— A) in an interval, with zero exterior condition, for 
complete Bernstein functions ip such that p'lpff) converges to infinity as £ —> oo. This 
extends previous results obtained by the authors for the fractional Laplace operator 
(ip (£) = p a>/ ' 2 ) and for the Klein-Gordon square root operator (ip(£) = (1 + £ 2 )F 2 — 1). 
The formula for the eigenvalues in (—a, a) is of the form A„ = ip(pf) + 0(pf), where p n is 
the solution of fi n = — -’d(ii n ), and r d(p) £ [0, -|) is given as an integral involving ip. 


1. Introduction and statement of the results 


This is the final one in the series of articles where asymptotic formulae for eigenvalues 
of certain pseudo-differential operators in the interval are studied. The fractional Laplace 
operator (—A) a//2 was considered in [18] for a = 1 and in [20] for general a G (0, 2), while 
in [IB] the case of the Klein-Gordon square-root operator (—A + l) * 1 / 2 — 1 was solved 
(A dentotes the second derivative operator, the Laplace operator in dimension one). 
In the present article we extend the above results to operators ip{— A), where ip is an 
arbitrary complete Bernstein function such that £ip'(£) converges to infinity as £ —> oo. 

Let A n denote the nondecreasing sequence of eigenvalues of ip(— A) in an interval D = 
(—a, a), with zero condition in the complement of D. Furthermore, for p > 0 define 
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We note that G [0, |) and = 0(A) as /i —* oo. Finally, let p n be a solution of 
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We remark that the solution is unique for n large enough, and 

~ 2^ ~~ a^( n7r )/(2a) + O(-). 

The following is the main result of the present article. 

Theorem 1.1. If ip is a complete Bernstein function and lim^oo £ip'(£) = oo , then 


= rp{f&) + 0(A) 


as n —» oo. 


( 3 ) 


In many cases, p n can be approximated with more explicit expressions, at the price of 
a weaker estimate of the error term. We provide two examples. 

Example 1.2. Let ip{ff) = £ a//2 + 2 , where 0 < f3 < a < 2. Then (see Example 12.111 ) 
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Example 1.3. If 0 is regularly varying at infinity with index j £ (0,1], then one has 
lim^oo (see (fT51) ). Therefore, 


l^n 


nix 

2a 


(2—a)ix 
8 a 


+ 0 ( 1 ), 


and, using Karamata’s monotone density theorem, one easily finds that 




Remark 1.4. The moderate growth condition lim^oo £0'(£) = oo is satisfied by all 
regularly varying functions with positive index. It is not satisfied by a slowly varying 
complete Bernstein function 0(£) = l°g(l + £)■ There are, however, slowly varying 
functions which do satisfy the moderate growth condition, for example, 

,/rt f°° Z log zdz 

m = l ' 

which is asymptotically equal to |(log£) 2 as £ —* oo. 

Remark 1.5. Numerical simulations for 0(£) = £“/ 2 (using the results of [IT]) strongly 
suggest that the error in (J3]) is in fact of order 0{\). There is no numerical evidence 
for more general functions 0. We believe that at least when 0(£) = y/£, one can use a 
method applied in a somewhat similar problem in Dm to obtain a version of (|3]) with an 
additional term in the asymptotic expansion and an improved bound of the error term, 
but this is far beyond the scope of the present article. 


We point out that relatively little is known about \ n . Most results, including all listed 
below, cover also higher-dimensionsinal domains, but provide significantly less detailed 
information. Extension of Theorem o for higher-dimensional domains seems out of 
reach with the present methods. 

Best known estimates of A n , proved in |8|, are given in terms of the corresponding 
eigenvalues A0 of the Laplace operator —A, namely 

CWA?) < A„ < V(A?); 


a more direct statement for the case of an interval is given in (TT21) below. First term 
of the asymptotic expansion of A n , namely A n ~ 0(A^), is given in many cases in 0. 
This result follows by a Tauberian theorem from the asymptotic expression for the trace 
e~ tXn as t —» 0 + . 

Second term of the asymptotic expansion of the trace has been found in [2] [3j for 
(— A)“/ 2 , in [H 2?>! for (—A + 1) Q//2 — 1, and finally in [7] for a rather general class of 
isotropic Levy processes with unimodal Levy measure, satisfying some mild regularity 
conditions. Tauberian theory is, however, insufficient to obtain a result similar to Theo¬ 
rem 11.11 from the two-term expansion of the trace. To the knowledge of the authors, no 
results of this kind are known for domains other than intervals, with the only exception of 
the well-studied classical situation of the Laplace operator —A. The only related result, 
proved in 1121 . provides a two-term asymptotic expansion of Cesaro means jj ^2^=1 ^ n ^ or 
(—A )"/ 2 using the methods of semi-classical analysis. 

The proof of Theorem 1 1.1 1 is based on the explicit expression for the generalised eigen¬ 
functions of the operator 0(—A) in the half-line, found in |18| for (—A) 1 / 2 , and in [TUI : Z2] 
for 0(—A) for a general complete Bernstein function 0. The asymptotic expression (J3J) 
for (—A) Q / 2 simplifies to 


a» = (g - +00 


8 a 


as n —y oo, 
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because d /t = ( 2 ~ g a ) 7r . As mentioned above, this was proved for a = 1 in [T8], with constant 
1 in the asymptotic notation O(-), and for general a 6 (0,2) in [20 j , with a rather big 
constant in the term 0(1). A very careful estimate of [T6] yielded a version of ([3]) uniform 
in a > 0 for the operator (—A + 1) 1//2 — 1. In the present article we do not pay attention 
to the constant in the asymptotic term O(-). All our estimates are, however, explicit, 
and so it is theoretically possible to trace the dependence of this constant on a and 0. 

To facilitate the reading of the article, we sketch the main idea of the proof. The 
generalised eigenfunction of 0(—A) in the half-line (0, oo) corresponding to the eigenvalue 
0(/i 2 ) is given by an explicit formula F /t (x) = sin(/ix + 10) — G^x), where G M is the 
Laplace transform of a certain non-negative measure (here ‘generalised’ essentially means 
‘not square integrable’). We construct approximation (p n to eigenfunctions of 0(—A) in 
(—a, a) by interpolating between Tj t (a + x ) near —a and ±F fl (a — x) near a. In order 
that the sine terms agree, we need to set /i = /i n defined in ([2]). Due to non-locality 
of 0(—A), (p n is not an eigenfunction; nevertheless, we show that the L 2 (D) distance of 
0(—A)0 n and /i n Cp n does not exceed O(-) (Lemma 13.611 . This is sufficient to prove that 
there is some eigenvalue Xk( n ) within the O(-) range from 0(/x 2 ). Using the assumption 
that 00'(£) diverges to inhnity as £ —» oo, one easily hnds that the numbers k{n ) are 
distinct for sufficiently large n. It remains to estimate the number of eigenvalues A j not 
counted as Xk( n )'- this turns out to follow from an estimate of the trace (Lemma 13.11H . 

We conjecture that ([3]) holds for arbitrary complete Bernstein functions, without the 
moderate growth condition lim^oo £0'(£) = oo. Note that, however, if this growth condi¬ 
tion is not satisfied (for example, when 0(£) = log(l + 0) and a is large enough, then one 
cannot expect that the numbers k(ri) are distinct. Therefore, an extension of Theorem ll.il 
to general complete Bernstein function would require a completely different approach. It 
is also natural to expect that (]3j) holds for more general functions 0, for example, for all 
Bernstein functions 0 satisfying the growth condition. However, no expressions for the 
generalised eigenfunctions are known unless 0 is a complete Bernstein function, and 
so our approach cannot currently be used in this case. 

The method described above has been designed in [TSj and sucessfully used in [20] 
and [TB]. The core of the argument remains the same in the present article. Nevertheless, 
proving Theorem 1 1.1 1 in this generality requires rather non-obvious estimates of 
and G^x), as well as many other modfications; for example, the trace estimate in the 
final part of the proof needed some improvements. 

The remaining part of the article is divided into two sections. In Preliminaries, we recall 
definitions and basic properties of complete Bernstein functions (Section 12.Hi and the op¬ 
erator 0(—A) in full space R and in (bounded or unbounded) intervals (Sections 12.2H2.4I) . 
We also recall known properties of the eigenvalues X n (Section 12.51) and the generalised 
eigenfunctions F^x) (Section 12.61) . Finally, we prove the necessary estimates of (Sec¬ 
tion [2ZZJ) and G0(0) (Section l2.8p . The proof of Theorem 1 1.1 1 is given in Sectional which is 
divided into five parts. Pointwise estimates for 0(—A) (Section 13. lj) are taken from [IB] . 
Construction of approximations to eigenfunctions (p n is followed by a technical lemma, 
which asserts that 0 n is in the domain of 0(—A) in the interval (Section [372]) . An estimate 
for 0(—A)0 n (Section 13.31) follow then easily from the results given in Preliminaries. This 
is used to find estimates for the eigenvalues Xk( n ) (Section 13.41) . We conclude the proof 
by showing that k{n ) = n for n large enough (Section 13.51) . As a side-result, we obtain 
some properties of the eigenfunctions, listed in the final Section 13.61 
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2. Preliminaries 


All functions considered below are Borel measurable. For p £ [1, oo) and an open 
set D C R, the Lebesgue space L P (D ) is the set of functions / on D such that 
||/||lp(u) = (Jo \ f(x)\ p dx) l ! p is finite, and / £ L°°(D ) if and only if the essential supre- 
mum ||/||l oo (d) of |/(x)| over x E D is finite. The space of smooth functions with compact 
support contained in D is denoted by Cf°(D). By C 0 (D) we denote the space of con¬ 
tinuous functions in R which are equal to 0 in R \ D and which satisfy the condition 
lim x ^ ±00 /(x) = 0 . 

The Fourier transform of a function / £ L 2 (R) is denoted by T f. If / £ L 2 (R)flL 1 (R), 
then J-’f(C) = f(x)e~^ x dx. The Laplace transform of a function / is denoted by £f, 
£/(£) = f(x)e~^ x dx. Symbols x, y, z are used for spatial variables, while £, 77 , // 
typically correspond to ‘Fourier space’ variables. 

We sometimes use standard asymptotic notation: we write f{n) = 0(g(n)) if 
limsupn^, \f(n)/g(n)\ < 00 , and f(n) = o(g(n)) if lim^^ \f(n)/g(n)\ = 0 . 


2.1. Complete Bernstein functions. In this section we recall several classical defini¬ 
tions. A function f(x) on ( 0 , 00 ) is said to be completely monotone if (— l) n f^ n \x) > 0 
for all x > 0 and n — 0,1,2,... By Bernstein’s theorem ( |25l Theorem 1.4]), / is com¬ 
pletely monotone if and only if it is the Laplace transform of a (possibly infinite) Radon 
measure on [0, 00 ). If / is nonnegative on (0, 00 ) and f is completely monotone, then / 
is said to be a Bernstein function. By Bernstein’s theorem, Bernstein functions have the 
representation 

f(x) = cx + c + f (1 — e~ zx )M(dz) (4) 

J (0,oo) 


for some c, c > 0 and a Radon measure M such that ^ min(z, l)M(dz) < 00 . The 
above formula extends to complex x such that Rex > 0 , and defines a continuous function 
holomorphic in the region Rex > 0. 

If the measure M in (J4j) is absolutely continuous with respect to the Lebesgue measure, 
and the density function is completely monotone, then / is said to be a complete Bernstein 
function. One easily verifies that in this case 


f( x ) 


If x m(dz 
cx + c- 1 — /- 

7T 7(0, oo) X + Z Z 


(5) 


for some c, c > 0 and a Radon measure m such that /( 0 ^ min(l/z, l/z 2 )m(dz) < 00 . 
The above formula defines a holomorphic extension of / in the region C \ (— 00 , 0]. 

Bernstein and complete Bernstein functions appear in a number of different areas of 
mathematics. For more information on these objects, we refer the reader to [25]. 

We will need the following technical result, proved in the Appendix. 


Lemma 2.1. Let f is a complete Bernstein function with representation (15|) . Let g be a 
holomorphic function in {w E C : | Argw;| < C\\ (with 0 < C\ < such that g(x) is 
real for x > 0, and let h be a continuous function on (0, 00 ). Denote 

G(y)= sup \g{z)\, H(y) = sup \h(x)\, 

i// 4<|2:|<4?/ 2//4<:r<4y 

|Argz|<Ci 

and suppose that 

POO 

G{x)H{x) < C 2 min(x _1 , x -2 ), C 3 = (1 + y)G(y)H(y)dy < 00 

J 0 
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Im(/(—e l£ x)g(e l£ x))h(x)dx 


( 6 ) 


(7) 


for x > 0. Then 

r roc 

/ g(x)h(x)m(dx) = lim — 

J (0,oo) e^0+ 7T _ ^ 

= lim — / Im(/(—e l£ x))g(x)h(x)dx. 

s^o+ 7r J 0 

Following [EU22], if ^ is a non-constant complete Bernstein function such that if( 0) = 
0 , and g > 0 , we denote 

/ - 1 ~ £/h 2 

i -m/w) 

for (6 C \ ((—oo, 0] U {g 2 }), and if^g 2 ) = if(g 2 )/(g 2 if'(g 2 )). We also let 


^ f (0 = exp ( - 




7T J o £ 2 + r 


log^(r 2 )dr 


for £ G C with Re£ > 0. Then if^ is a complete Bernstein function, if' extends to 
a complete Bernstein function, and we have the Wiener-Hopf identity if' (£)if' (—£) = 
t^(-e 2 ) for (6C\R, see, for example, |T9i, Proposition 2.19 and Lemma 3.8]. Finally, 
we denote if^ = (if^Y- 

In principle we could extend the definition of if^ to general non-constant complete 
Bernstein functions if, so that = (1 — £ 2 //i 2 )/( 1 — (if(f) — if(0))/(if(g 2 ) — if( 0))). 

All results proved below hold true with this definition. However, to keep the notation 
simpler, we will typically assume that ^>(0) = 0. For brevity, we also denote if(oo) = 

liin^ 

OO V»(0 e [0,oo]. 

Unless otherwise stated, in what follows we assume that if(£) is a non-constant complete 
Bernstein function which satishes ^(0) = 0, that is, c = 0 in representation ([5]) for if. 


2.2. Analytical definition of the operator if (—A) in R and in intervals. In this 
section it is enough to assume that if is an increasing, nonnegative function on [ 0 , oo), 
which satisfies 


l + if{f + i])<C(l + r]) a (l + if(f)) ( 8 ) 

for all £ > r) > 0 and some C, a > 1. When if is a complete Bernstein function, then (j5|) 
holds with a = l and (7 = 1, because if(f + rj) < if(f) + if{rj) < if(f) + (7(1 + rj). 

The operator A = if (—A) is an unbounded, non-local, self-adjoint operator on L 2 (R), 
defined as follows. The domain T>{A) of A consists of functions / G L 2 (R) such that 
(1 +if(f 2 ))J r f(£,) is square integrable. Clearly, V(A) contains C£°(R). For / G V(A), 

FAf(o =v^ 2 )^m 

In other words, A is a Fourier multiplier with symbol if((, 2 ). This explains the notation 
A = if (—A): the second derivative operator A is a Fourier multiplier with symbol — £ 2 . 
Furthermore, by Plancherel’s theorem, A is positive-definite. 

Let T>{£) denote the space of / G L 2 (R) such that (1 + '0(£ 2 )) 1 // 2 7 r /(£) is square 
integrable. For /, g G T>(£) the quadratic form £(f,g) associated to A is defined by 

£U,g) = TJj, 

The inner product £\ (/, g) = (/, g)+£(f, g) makes T>(£) into a Hilbert space. If / G F>(A), 
then £(f,g ) = ( Af,g ), and T>(A) is a dense subset of the Hilbert space T>{£). 

Let D be an open subset of R. The following definition states that the operator Ap is 
the Friedrichs extension (or the minimal self-adjoint extension ) of the restriction of A to 
C?{D). 
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Definition 2.2. The domain D(£ d ) of the form £ D is the closure of Cf°(D) in the Hilbert 
space T>(£), and £n(,f\ g) = S(f,g) for g G V(£ D ). The operator A D is associated 
to the form £q: f G T>(£ D ) is in the domain V(A^) of Ap if and only if there is a 
function A^f G L 2 (D ) such that £(f,g) = (Ar>f,g) for g G T>(£d) (or, equivalently, for 
geC?(D)). 

The following result is well-known in the context of general Dirichlet forms and gen¬ 
erators of Levy processes, see [131 [2i1j for more general results in this direction. For 
completeness, we provide a short proof. 

Proposition 2.3 (see [XBj Proposition 2.2]). If D is a bounded interval, then f G T>(£o) 
if and only if f G T>(£) and f = 0 almost everywhere in R \ D. 

Proof. By definition, if / G T>(£d), then / G T>(£) and / = 0 almost everywhere in 
R \D. Let / G T>(£) and / = 0 almost everywhere in R \D. The result follows from the 
following claim: there is a sequence f n G Cf°(D) such that 

Si(f n -fjn-f) = ^J 

converges to 0 as n —* oo. 

Let h n G CP°(R D ) be an approximation to the identity such that h n (x) = nh(nx ), 
h(x) > 0, f R h(x)dx = 1 and h(x) = 0 for x ^ (—1,1). Note that h n is zero outside 
(--,-). Let 

g n (x) = h n * f(x ), f n (x) = g n ((x - b n )/a n ), 

where (x — b n )/a n maps the ^-neighbourhood of / into /, with a n > 1 , lirn^oo a n = 1 
and lim .,,^00 b n = 0. Observe that f n G Cf°(D) and 

FfniO = a n e~ ibn ^J r g n (a n tf) = a n e- ib ^h(S(a n 0) -F/M- 

Since /,g G L^R), Tf and Th are continuous. Furthermore, ThifS) = 1 and |J-7i(£)| < 1 
for £ G R. It follows that Tf n converges pointwise to Tf , and for n large enough, 

|^/n(0l<2|^/(fln0l 

for all (gR. Hence, if u(f) = (1 + \J- f(f) | 2 , then for n large enough, 

(i+^ 2 ))rao - i 2 < 2(1++i^/(oi 2 ) 

< 4 u(a n f) + 2 u(f) 

for all f. By the assumption, m(£) is integrable. Therefore, the family of functions 
f n (t;)—J 7 f (£,)\ 2 is tight and uniformly integrable. By the Vitali’s convergence 
theorem, £\ (f n — f, f n — f ) converges to 0 as n —* oo, as desired. □ 

We remark that the above result in general fails to be true for arbitrary open sets D. 
It is, in particular, not true when D = R \ {0} and f>(f) = f a ^ 2 with a G (1, 2], 

2.3. Markov semigroup generated by A. From now on, if is a complete Bernstein 
function. The operator —A = —-0(—A) generates a strongly continuous semigroup of 
self-adjoint contractions 


T{t ) = exp(— tA), 
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where t > 0. Note that T(0) is the identity operator, T(t) is the Fourier multiplier with 
symbol exp(—h0(£ 2 ))j and 


V(A) = <j / £ L 2 (R) : lim 


2 , 0 ,. r(^)/ — / _ r2 , 


exists in L (R) 


—Af = hm 
t->o+ 


T{t)f - f 


(9) 


For t > 0, the operator T(t) is a convolution operator: for all / £ L 2 (R), 





y)T(t; dy), 


( 10 ) 


where T(t]dx ) is a sub-probability measure with characteristic function exp(—H/;(£ 2 )) 
and total mass . Furthermore, 

T(t- dx) = e- t ^ oo) 5 0 (dx) + T(t; x)dx, 


where T(t\x) = T(t ; —x) is a decreasing function of x > 0 (see [23]). Hence, T(f) is a 
Markov operator , and formula (fTUl) defines a contraction on every L P (R) ( p £ [l,oo]), 
and also on C*o(R). In each of these Banach spaces, the generator of the semigroup T(t) 
is dehned in a similar way as in (|9]1; for example, 


V(A-,C 0 (R)) 


-Af 


f e C 0 (R 
T(t)f - f 


hm 

t-> o+ 


T(t)f - f 


exists in C'o(R) 


lim 

t->o+ 


t 


Since the above definitions of Af are consistent on the intersections of domains with limits 
in different function spaces: L P (R) for p £ [l,oo] or C'o(R), we abuse the notation and 
use the same symbol —A for the generator of the semigroup T(t) in any of these spaces. 
Observe that (^(R) is contained in V(A,L P ( R)) ( p £ [1, oo]) and in D(H; C 0 (R)), and 
it is the core of A in each of these Banach spaces except L°°(R) (see Bill). Whenever 
we write T>(A), we mean V(A] L 2 (R)). 

If 'f(^) has the representation given in ((3]), then for / £ C^°(R) we have 

/ OO 

(f(x) - f(y))u(x-y)dz, (11) 

-OO 


where by the subordination formula, 


u(z) 


J_ /‘ - k | C i/ 2 m(dC) 

2vr./(o,oo) C 1/2 


and ‘pvj’ denotes the Cauchy principal value integral: 


P v / (/0) ~f(x + z))v{z)dz 


lim / (/(x) - /(x + z))v(z)dz\ 

£ ^ 0+ J R\(-e,e) 


see, for example, m Proposition 2.14], 


2.4. Markov semigroup generated by Ad- Let D be a (possibly unbounded) interval. 
The operator — Ad generates a strongly continuous semigroup of operators 

T D {t) = exp(—C4 d). 

The operators To(t) are given by 

T D {t)f{x) — [ f(y)T D (t-,x,dy ), 

J D 
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where 

T D (t; x , dy) = e~ t ^ oo) 8 x (dy) + T D (t ; x, y)dy. 

It is known that 0 < T D (t]x,y) < T(t,x — y), and we let T D (t',x,y) = 0 whenever 
x (fc D or y £ D. Hence, Tjj[t) form a contraction semigroup on each of the spaces L P (D) 
(p G [1, oo]), and if ip is unbounded, then also on C 0 (D ) (see |S1 HH1123|)• The generator of 
each of these semigroups is again denoted by — An , and it acts on an appropriate domain 
V(A D -L p ) or V(A d -Co). 

2.5. Spectral theory for Suppose that D is a bounded interval and that 

exp(— 2tip{tf 2 )) is integrable for some t > 0. Then To(t;x,y) is a Hilbert-Schmidt kernel, 
and so To(t) is a compact operator on L 2 (D). Hence, there is a complete orthonormal set 
of eigenfunctions ip n G L 2 (D) of Tr>(t). By strong continuity and the semigroup property, 
the eigenfunctions do not depend on t > 0 , and the corresponding eigenvalues have the 
form e~ tXn for all t > 0 , where the sequence X n is nondecreasing and converges to oo. 

By translation invariance, with no loss of generality we may assume that D = (—a, a). 
By symmetry, T D (t',x,y) = T D {t ; —x, —y), and hence the spaces of odd and even L 2 (D) 
functions are invariant under the action of T D (t). Therefore, we may assume that every 
(p n is either an odd or an even function. The ground state eigenvalue X\ is positive 
and simple (unless ip is constant), and the corresponding ground state eigenfunction has 
constant sign in D ; we choose it to be positive in D. The functions ip n are also the 
eigenfunctions of Ad (because —Ad is the generator of the semigroup To(t)), and X n are 
the corresponding eigenvalues. 

No closed-form expression for X n and ip n is available, except when ip(f) — cf + c. By 
a general result of [5] (see Theorem 2.3 therein), X n ~ ^((^y) 2 ) as n —> oo (the original 
statement includes only the case when ip(£) ~ for some a G (0,2), but it can be 
easily extended to more general ip). Best known general estimates of X n are found in [HJ 
Theorem 4.4], where it is proved that: 



Note that the upper bound in (1X21) follows relatively easily from the operator monotonicity 
of ip: the form associated to Ad is bounded above by the form of iP(-Ad), and the 
eigenvalues of the latter are equal to ^((ff) 2 )- The proof of the lower bound is more 
intricate. 


2.6. Spectral theory for H( 0jO o)- The spectrum of Ad for an unbounded interval D is 
continuous. When D — R, then Ad = A takes diagonal form in the Fourier space, and e l ^ x 
(£ G R) are the L°° eigenfunctions of A. Similar eigenfunction expansion was obtained 
for the half-line using an appropriate version of the Wiener-Hopf method in | K[[ [22]. Due 
to translation invariance and symmetry, it suffices to consider D = (0, oo). 


Definition 2.4. Suppose that ip is a non-constant complete Bernstein function such that 
-0(0) = 0. For x, n > 0, let 


F^x) = sm(px + - G m (x), 

where G [0, |) and G^ is a completely monotone function on (0, oo). More precisely, 




d 


log 


ip'(lJ, 2 )(p 2 — r 2 


r 2 — y 2 ip{d 2 ) ~ ip(r 2 ) 


■ dr 


7r 
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(as in ©), and G M is the Laplace transform of a measure 7 ^, 


G„(x) = C^^x) = 




' (0,oo) 


7 MO, 


with 


= — lim Im 

7r e-> 0 + 


f 




x exp 


x 


K 


£ 


f 2 + r 2 


g - ^(r 2 ) 1 4 


for p,f,x > 0 . 

Equivalently, F^(x) is dehned by its Laplace transform: for (eC with Re£ > 0, 


/*CO 

£F m (0 = / = 

7o 

see |22l Theorem 1.3] and J 

= 




h 2 + f 


exp — 

7T 


£ 


log W)hlxL) rfr 


£ 2 + r 2 ^(/i 2 ) ~~ ip(f 2 / 

Theorem 1.1]. We have the short-hand expressions 

h V>J (0 


vwr 

^m = Arg 


7m(^ 0 = - lin \ 


A* Im^(-e i£ £ 2 ) 




vr^o+^ + e 2 ^(/r 2 )^(0 

again see H3 Remark 4.12] (note the typo in formula (4.14) therein) and [221 formu¬ 
lae (2.17)-(2.19)]. The expressions for j^df) given above are slightly different than 
in mm, so we provide a short justification. By Lemma [2. II and the identity 

p cos tdfj, + £ sin 


£G *(0 = 


h 2 +£ 2 


^( 0 , 


we have 


7/x(d0 = 7 limIm(£G #1 (e *®0)df = ~ lim Im(£F M (e lE £))d£. 

£— ^0 ' 

The expression for £F^(£) and the Wiener-Hopf identity (^)V’j t ( — 0 = ^(—■£ 2 ) 


give 


7m (df) = - lim Im I 
7 r e-s>o+ 


= — lim 


£ 


*1’n 


~—2ie c 2 


a 


M 2 + e 2ie £2 t 

Im'tpJ— e~* £ £ 2 ) 


df 




7T E^0+ /l 2 + £ 2 ^m(/^ 2 )^m(0 

as desired; here we used Lemma m again. 

We extend the dehnition of F /; and G fl to R so that F^ix) = G lt (x) = 0 for x < 0. The 
functions F /; (// > 0) are L°° eigenfunctions of and play a similar role for Ad as the 
Fourier kernel e l ^ x (£ G R) for A. This is formally stated in the following result. 


Theorem 2.5 ([T9J Theorem 1.1] and [(221 Theorem 1.3]). The functions F M are L°° 
eigenfunctions of A^ 0,00)?' the corresponding eigenvalues are if{p 2 ). The operator A( 0)OO ) 
takes a diagonal form under the integral transform with kernel F tl {x) . More precisely, let 


n f{iT) = / f(x)F fl (x)dx 
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for f G L 2 ((0, oo)) D L 1 ((0,oo)). Then (-) 1 ' /2 II extends to a unitary mapping on 
L 2 (( 0, oo)), such that for f G L 2 ((0, oo)), 

« (l + ^(/i 2 ))n/(/i) g l 2 ((o, oo)), 

and if f G X?(AL(o,oo)) 7 i/ien 

n(Ao,oo)/)M = W)n/0A n(r D (t)j) = 

In this article we only use the first part of the above result, namely, that Ff are the 
L°°((0, oo)) eigenfunctions of A( 0 ,oo)- We remark that a similar eigenfunction expansion is 
available for D = R.\ {0}, see dSlEU, and there are no other known explicit expressions 
for the eigenfunctions of Ad unless D = R or if(f) = cf + c. 


2.7. Estimates of Recall that according to HD. Definition 12.41 and [ 191 Proposi¬ 
tion 4.16], 


= 


T 


T 


log 


1-z 2 


log 


y>y)(h 2 -* 2 

lf(p 2 ) — lf(p 2 Z 2 


ds 


dz. 


(13) 

(14) 


z 2 {if(p 2 /z 2 ) — if(p 2 )) 

We remark that if if is regularly varying at infinity with index a G (0,2], then, by 
dominated convergence, 

1 /•! 1 


lim TP = — lim 

[i —^oo 7T fi—> oo 


l 

7T 


~1 


log 


log 


1 - lf(p 2 Z 2 )/lf(p 2 ) 

z 2 (if(p 2 /z 2 )/if(p 2 ) — 1 ) 


dz 


l~z c 


dz = 


2 — a f 1 — log z 


1 — z 2 


dz = 


(2 — o)7T 


(15) 


I o 1 — Z 2 z 2 (z~ a — 1) 7T Jo 

see [ITT . Example 6.1]. By [191 Proposition 4.17], dominated convergence can be used to 
differentiate the right-hand side of (fT3l) in p > 0 under the integral sign. This yields 

1 'p 2 if\p 2 )-p 2 z 2 if'{p 2 z 2 ) {p 2 /z 2 )if\p 2 /z 2 )-p 2 if'{p 2 )' 


y = _ 2 _ , x 
dp np J 0 1 — z 2 


if(p 2 ) — if{p 2 z 2 ) 


if(p 2 /z 2 ) — if{p 2 ) 


dz 

(16) 


for all p > 0. In this section we prove two properties of d^ that are needed in the 
remaining part of the article. First, we find estimates of d^ that imply that the lower 
limits of as p —» 0 + or p —* oo do not exceed ^ (Lemma 12.91) . Next, a simple estimate 
of j^d/j, is found (Lemma 12.10]) . 

By [22, Proposition 4.3], we have the following general estimate of d^: 



-W(0 \ 

no ) 




sup 

«>0 




7r 

4 


for all p > 0. Furthermore, the supremum in the upper bound is always not greater 
than 2. If if is a Thorin-Bernstein function (see [25]), then one easily checks that the 
supremum is in fact not greater than 1, and therefore < j. Below we find more refined 
bounds for d^. 

By [22, Proposition 4.5], 




7T 




arcsm 


\[Q + 


arcsm 


Q 


l -p 


arcsm" 


PQ 

l-P 


< d„ < - 
- 11 ~ 2 


arcsin yfP (17) 
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with 

= /iW) = -hW) 

2 ) ’ WW) 

(note that the factor - is missing in the lower bound in the original statement). By the 
same argument as in the proof of the lower bound of m Proposition 4.5] (using the lower 
bound for Vt\(A 2 £ 2 ) and the upper bound for ^(A 2 /£ 2 )), one easily shows that, with the 
same P and Q, 

^ — — ( arcsin 2 y/l — Q + arcsin 2 \ 1 — ^ — arcsin 2 \ 1 — ] . (18) 

4 7r \ VI — P V 1 — P J 

One can also verify that this bound is always at least as good as the upper bound of (fT7j) , 
with equality when P + Q = 1. 

The following technical result states that P+Q < 1. This in fact follows indirectly from 
the proof of |22l Proposition 4.5] (note that the right-hand side of (1T5]) is not well-defined 
when P + Q > 1), but we choose to give a simple, direct argument. 


Lemma 2.6. If if is a non-constant complete Bernstein function, then 

V>'(£) “ V’(f) 

Proof. The lemma is equivalent to the inequality 

-mono < w(om) - mo)- 


Assuming if has the representation (]5|) . we need to prove 


£ c£ + c + - 
7r 


< c + 


f £ K ds ) 

qo,oo) £ + s s 


p,(ds) 


n(ds) 


(0,oo) (£ + S Y 

i r 

c + 

r (0,oo) 


£ 2 h(ds) 


(£ + *) s 


* 7(0,oo) (£ + S ) 2 * / V * 

This follows by simple integration from the following bounds: 0 < cc, 

S - £ 2 


£(c£ + c) 


(£+«) 


< c 


3 — 


+ c 


and 

£ 


£ 


«2 


+ 


£ 


Sl 


< 


(£ + s ) 2 ' " (£ + s ) 2 ’ 

£ 2 


Sl 


+ 


s 2 


£ 2 


£ + Si (£ + s 2 ) 3 £ + s 2 (£ + Si) 3 / (£ + Si ) 2 (£ + s 2 ) 2 (£ + s 2 ) 2 (£ + Si ) 2 ’ 

the last two inequalities are easily proved by direct calculations. □ 


Lemma 2.7. The left-hand side of (H71) is decreasing in P G [0,1 — Q]. The right-hand 
side of (1151) is increasing in P G [0,1 — Q\. 


Proof. Let P = 1 — -42_ = -A— 5 s£ [0,1 — Q]. Note that P increases with increasing .s, 
and the left-hand side of (1T71) is equal to 

^(arcsin 2 \fQ + arcsin 2 \J s + Q — arcsin 2 Vs). 

Since arcsin 2 y/s is convex, the above expression is increasing in s. In a similar way, with 
P = 1 — s E [0,1 — Q\, the right-hand side of (H5|) is equal to 

j — ^(arcsin 2 a/1 — Q + arcsin 2 \Js — arcsin 2 V s + Q)i 
which is again an increasing function of s, but now P decreases with increasing s. □ 
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Substituting P = 0, we obtain immediately the following elegant result. 


Corollary 2.8. If if is a non-constant complete Bernstein function such that if(0) = 0, 
then 


arcsm 


-p 2 if"(p 2 


7T 


2 if'(p 2 


< ^ < - 
_ n 2 


arcsm 


'1 + 


7T 


p 2 lf"(p 2 

WW) 


Lemma 2.9. If if is a non-constant complete Bernstein function such that if( 0) = 0, 
then 


If if is unbounded, then also 


3 7T 

lim inf < — . 

/x—> o+ 8 


3yr 

lim inf d u < — . 


Proof. Suppose that lim inf /i ^ 0 + Then there are po > 0 and q G (0,1) such that 

> | — ^ for p G (0, /iq)- By Corollary 12.81 


arcsm" 


1 + 


p 2 if"{p 2 ) ^ qn 

w 


21 ^ 


16 


for p G (0, /U 0 ) 5 and hence 


W'V) 

if'(p 2 ) 


>2 — 2 sin 


. 


for p G (0,/io)- If « denotes the right-hand side, then a > 1. By integration (see fl5l 
Lemma 2.2]), we have if\p 2 )/if\pfy > {pi/p 2 ) a for all p G (0,/i 0 ), which contradicts 
integrability of if' at 0. This proves the first statement of the lemma. 

In a similar manner, if lim infthen there are /i 0 > 0 and q G (0,1) such 
that | ^ for p G (/io, oo). Again this implies 

-/< W) >2 -2 (sin " 


lf'{p 2 ) 

for p G ( po,oo ). If a denotes the right-hand side, then a > 1, and by integration, 
if'(p 2 )/if'(pi) < (pl/p 2 ) a for all p G (/io,oo). This implies integrability of if' at oo. □ 

We conjecture that the above lemma holds with ^ replaced with An example of a 
complete Bernstein function if for which the set of partial limits of d^ as p —> 0 + is equal 
to [0, f] is given in [221 Section 7.5]. 

Lemma 2.10. If if is a non-constant complete Bernstein function such that if( 0) = 0, 
then for all p > 0 , 


dd„ 


dp 


3 

< - . 

P 


Proof. By (TT 6 l) and the Cauchy’s mean value theorem, for some f z G (p 2 z 2 ,p 2 ) and 
ii/z e (h 2 ,h 2 /^ 2 ) (where z G (0,1)), 


dd h 


dp up J 0 1 


np Jo i - - 2 v y(y) 


wit *)+y( 6 ) Zi/znti/z )+y 

t*r{tz) ti/sV'(ti/z) 


tPitl/z 


dz 


y(6 A 


dz. 
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By (J5]), 0 < — CV ;// (0 < 2'0'(£) and 0 < £ 2 /M 3 )(£) < 6-0'(£). Hence, 

P d (mt)\ e^ 3 \o -mo r r fi fil 

^(o ^(o moy [ ’ J ' 

Furthermore, ~ £i/z'4’"(£i/z)/^'(t,i/z) 6 [-2,2]. It follows that 


dW 


d/i 


< — 


7r/i Jg 1 — £ 2 


min 2 , 


£l/z 


d / ri/j"(r) 


— 


< 


71/1 J 0 1 — Z 

Recall that £i / 2 /£ 2 < z -4 . Hence, 

dtf 


min f 2 , 6 log ) dz. 


dr y 

dZ 

& 


dr dz 


d/i 

Since — log £ <2 — 1 , we have 

dd 


< 2 j' 1 min(2, — 24 log z) 


T/i Jo 


1-z 2 


d/i 


2 f 1 min(2,24(l/z- 1)) ^ 100 log5 - 48 log24 ^ 3 

8 / — d, — y 


do 


1-z 2 


□ 


7T/i 


d 


We conjecture that in fact — - < We close this section with the following 

simple example. 

Example 2.11. Let -0(£) = £"/ 2 + ^Z 2 , where 0 < f3 < a < 2. By a short calculation, 
//’(/i 2 ) — //;(/i 2 z 2 ) 1 ^1 — z“ ^ 1 \ /1 — z Q ^ z Q_/3 


-i 


z 2 (i/>(d7- 2 ) -V>(d 2 )) z 2 ~ a yl — z@ /i Q_ dy y 1 — z@ /i" _ d 
If we denote w = (1 — z a )/( 1 — z@), then 

- 0 (/i 2 ) — - 0 (/i 2 z 2 ) 1 n a ~^w + 1 1 f i _\_ 1 — £ q-/3 


z 2 (//;(/i 2 /z 2 ) — //>(d 2 )) z 2_Q /i" _/ % + z 0 —- 8 z 2_ “ y fi a -Pw + z“ _ dy 

As in the last equality of (1151) . we obtain 

q (2-a)7r 1 f 1 1 , / 1 - 2 ““^ \ , 

' ? "=— s — + ^ 7 0 iv^y 

Clearly, the integrand is nonnegative. Since log(l + s) < s, z a ~ 13 > z 2 and w > 1, 


1 -Z‘ 


log 1 + 


1 - z a ~P 

li a ~Pw + 1 


dz < 


< 


1 - z a ~P 


1 — Z 2 ll a ~Pw + Z a ~P 


dz 


1 1 - - 2 , 1 

■ dz = 


1 — z 2 fi a ~P n a ~P 


(2 — aW (2 — aW 1 

1 ’ <‘d^< 1 ; + 


7TyU' 


a— (3 


Therefore, 
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2.8. Estimates of F u( x )- In the remaining part of the article we will need the following 
simple estimate of CFp and a more refined estimate of G^. 

Lemma 2.12. If ip is a non-constant complete Bernstein function such that ip(0) = 0, 
then for all p > 0 and f such that Re£ >0, 


\jCF,(f)\<2V2 



n^ 2 -m 


Proof. Recall that (p 2 + f 2 )CF^(f) = p{ip^(p 2 )) lF ipjXf) is a complete Bernstein function 
of £, and hence by pT9], Proposition 2.21(c)] and [22] Corollary 5.1], 


l/i 2 + e 2 | \CF,(0\ < V2 (p 2 + |£| 2 )£F M (|£|) < 2v / 2p 


^(/i 2 ) -i/>( i^i 2 ) 


for all £ such that Re£ > 0. □ 

Lemma 2.13. If ip is a non-constant complete Bernstein function such that ip( 0) = 0, 
then for all p, x > 0 such that px ^ 1, 


1 iP(l/x 2 ) l iP'(pi) l-iP(p 2 )/{p 2 x 2 iP{l/x 2 )) 
11 X ~ TTX 1p(p 2 ) y -0(/l 2 ) 1 — V^(l /x 2 )/lp(p 2 ) 

In particular, if ip is unbounded, then 

lirnsup (pip(p 2 )GJx)) < . 

fi—> oo 


Proof. Recall that > ip^{ 0) = ip^(0) = 1. Hence, 


iMO < - / , , lim Im ip^—e l£ f 2 )df 

npy/ipuKP 2 ) ^ 0+ 


After a substitution f = yfs it follows that 


G,(x) = 


e < 


lim 


2np^/ip fl (p 2 ) e^o+ 


se 


-x^Ts Wyj-e i£ s)ds 


Since Xyfs e xF * < 2/(1 + x 2 s), we have 


G,(x) < 


-lim 

npxy/ip^p 2 ) £^o+ 


1 Im ipn{—e l£ s)ds ip^l/x 2 ) — 1 

l + x 2 s S ~ Trpxy/lp^p 2 ) ’ 


for the last inequality note that the integral converges to the integral term in the repre¬ 
sentation ([S]) for the complete Bernstein function ip and we have ip^( 0) = 1 (therefore 
the inequality becomes equality if ip^ contains no linear term, that is, if ip is unbounded). 
To prove the first statement, it remains to use the definition of ip^. The other statement 
of the lemma follows from the first one by the inequality £ip'(£) < ^(£). □ 


3. Proofs 

Throughout this section we implicitly assume that ip is a non-constant complete Bern¬ 
stein function such that ip(0) = 0, that is, c = 0 in the representation (J5J) for ip. By c 
and v we denote the constant and the measure in the representation (J5j) for ip. Finally, 
we let D = (—a, a) for some a > 0. 
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3.1. Pointwise estimates for the operator A. Recall that A = fi>(—d 2 /dx 2 ), and for 
/ G C£°(R) we have, as in (HIT) . 


Af(x) = -cf"(x ) + pv / (f(x) - f(y))u(x - y)dy 

J —oo 

poo 

= -cf"(x) + / (2 f(x) - f(x + z) - f(x - z))v(z)dz. 


(19) 


We denote the right-hand side by A fix') (with a calligraphic letter A) whenever the 
integral converges and, if c > 0, f" is well-defined. The following estimates of Af(x) 
are proved in [16] in the special case ip(£) = (£ + l) 1 / 2 — 1, but their proofs rely only on 
the symmetry, unimodality and positivity of the kernel function v. Note that in [16] the 
notation A o is used for A. 

Lemma 3.1 ([16j Proposition 4.1]). Let x G R, b > 0, and let g have an absolutely 
continuous derivative in [x — b, x + b). Then 


\Ag(x)\ < c\g"{x)\ 


sup \g"{y)\ 

y£(x-b,x-\-b) 


z 2 v(z)dz 


(l^(^)l + \g(y)\)v(y - x)dy. 


J TV\(x—b,x+b) 

As in usmaEa, for b > 0 we dehne an auxiliary function: 

0 for x G (—oo, — b ], 


q(x) = < 


(l/2){x/b + l) 2 for x G [—6, 0], 
1 — (l/2)(x/b — l) 2 for x G [0, b], 

1 for x G [b , oo). 


( 20 ) 


Note that q is C 1 , q' is absolutely continuous, 0 < q"(x) <l/b 2 (for x G R\{-6,0,6}), 
the distributional derivative q^ is a finite signed measure, and q(x) + q(—x) = 1. 

Lemma 3.2 ([IB) Proposition 4.2]). Let b > 0, let f G L 1 (R) ; and suppose that the 
second derivative f"(x) exists for x G [—b,b\ and it is continuous in [—6,6], Define 


M_! = 


I f{x)\dx, 


M 1 = sup | f{x) |, 


M 0 = sup | f (x) |, 

x(E[—b,b] 

M 2 = sup | f\x) |, 


xe[-b,b\ xe[-b,b\ 

Let q(x) be given by (j20|) . and define g(x) = q(x)f(x). For x G (—oo,0), we have 

|-4fi , (3 ; )| ^ C{b, '0)(il4_i + Mq + Mi + M 2 ). 

More precisely, for x G (—oo, —b] we have 


Mn 


c2b 


\Ag(x)\<f^ I z 2 v{z)dz + u(2b)M_i 


and for x G (— 6 , 0), 
\Ag(x)\ < M 2 c + 


M 0 2M\ 


z 2 uiz)dz + 2M 0 / v{z)dz + v(b)M_i. 
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3.2. Approximate eigenfunctions. Recall that D = (—a, a). Following [TB1 15] 120]. 
for n > 1, let Jx n be the largest solution of 


T $p,n 


nn 
2 ’ 


( 21 ) 


with dehned in (fljl ; this agrees with the dehnition of /i n in (|2]) , but we choose to use the 
notation ft n , so that all approximations are clearly distinguished from true values by the 
presence of a tilde. Although we are interested in large n only, note that by Lemma [2.91 
the equation a/i + ^ has a solution for all n > 1, and every such solution satisfies 


(n—l)7i 
2 a 


^ l~l j n ^ 


< 


niv 
2 a ' 


We remark that (12TT1 may fail to have a unique solution for n — 1 (for example, when 
a = 1 and -0(£) = £/(10 4 + £) + C/10 7 )- Nevertheless, if n > 3 and /i > then, 

by Lemma 12.101 


t ( a/i + > a -?> a -v >0 ’ 

and so the solution /!„ is in fact unique. 

We let 


An — V’(h'n)- 

In order to show that A n is close to some eigenvalue of Ad, we construct an approximate 
eigenfunction (p n of Ad, using the eigenfunctions F^ n {a — x), F^ n {a + x) for the one-sided 
problems corresponding to A(_ 00ia ) and A(_ a oo ). As in [T6l [151120] . we dehne 

<p n (x) = q(-x)F-^(a + x) - (-1 ) n q(x)F-^(a - x), (22) 

with the auxiliary function q dehned by (120T) . Here x G R, but we have (p n {x ) = 0 for 
x £ D, so that (p n is equal to zero in the complement of D. Clearly, (p n is continuously 
differentiable in D, (p' n is absolutely continuous in D, Cp" n exists in D \ {— b , b }, and (p" n is 
locally bounded in D. Note that A n depends on a and n, while (p n (x) depends also on 
b. We could fix b in order to optimise the constants (in many cases b — |a seems to be 
a reasonable choice), but since we do not track the exact value of the constants, we will 
simply indicate their dependence on b. Note also that (p n is not normed in L 2 (D ), its 
norm is approximately equal to \[7i (see Lemma 13.511 . 

The notation introduced above is kept throughout this section. 

The following result is intuitively clear, although its formal proof is rather long and 
technical. 

Lemma 3.3 (see [TB] Lemma 4.1]). We have <p n G V(A D ) and A D (p n (x) = Afi n (x) for 
almost all x G D. 


Proof. For brevity, in this proof we write fi = fi n and (p = (p n . The domain of Ad is 
described in Dehnition 12.21 we need to prove that (p G T>(£) and that (<p,Ag) = (Ap,g) 
for all g G We hrst verify the latter condition. 

Note that Ap(x) is well-defined for all x G D \ {—b, b}, since (p is smooth in D \{—b, b} 
and bounded on R. Let g G Since fJ is absolutely continuous in (—a, a), 

integration by parts gives 


-op"(x))g(x)dx = / <p(x)(—cg"(x))dx. 
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Furthermore, by the definition of A (see f|19p ). 

pa pa 

I A(p{x)g{x)dx — / (p{x)Ag{x)dx = 

(g(x + z)<p(x) + g(x — z)<p(x) — g(x)<p(x + z) — g[x)(p{x — z))v(z)dz ) dx. 


a / roc 


'—a \J 0 


We claim that the double integral exists. Then, by Fubini, it is equal to 0, and so 
( (p,Ag) = (A<f,g), as desired. 

Denote the integrand by I(x, z)u(z), and let £ — | dist(supp <7, R \D), so that supp g C 
(— a + 3e, a — 3e). When z > e, then \I(x,z)\ < 4||<^||x,°°(r) ||<?||l°°(r)- Suppose that 
0 G (0, e). If x ^ (—a + 2 e,a — 2 e), then I(x,z ) = 0. Otherwise, by first-order Taylor’s 
expansion of 7(x, z) around z — 0 (note that I(x, 0) = ^/(x, 0) = 0) with the remainder 
in the integral form, we obtain that 

f z d 2 

\I(x,z)\< J (z- s)—I(x,s)ds 

< ^ 2 (I|^|U° 0 (R)I|5 ,,, ||l°°(R) + ||||L°°((—a+e,a—e)) 115 1 11 L°° (R)) 


(recall that (p" is bounded in (—a + e, a — e)). We conclude that \I(x, z)u(z)\ < 
Ci((p, g) min(l, z 2 )u(z), which implies joint integrability of I(x, z)u(z). Our claim is 
proved. 

It remains to verify that <p G T>(£), that is, (1 + V’(£ 2 ))I^W’(0| 2 is integrable. Let 
f(x ) = q(a — x)Fp,(x), so that <p(x) = f(a + x) — (— 1 ) n f \a — x) (see (J22j)). It suffices to 
prove integrability of (1 + -0(£ 2 ))1 2 - 

Fix e > 0 and let q(x) = q(a — x)e £X . Since the distributional derivatives q, q' and 
q" are integrable functions, and the third distributional derivative of q(x) is a finite 
signed measure on R, the function q(x) has the same property. Therefore, J~q(£) and 
■Fq( 3 \£) = —i^Fq^) are bounded functions, and so |J r g(^)| < 0 2 (e, a, b)/(l + |£|) 3 . The 
Fourier transform of e~ £X Fp,(x) is equal to CF^e + *£), and the Fourier transform of 
f(x) = q(a — x)Fp,(x) = q(x)e~ £X F ^(x) is given by the convolution 

Ff(0 = ^:j J r q(£^s)CF fi (£ + is)ds. 


Suppose that £ > 0. To estimate |J r /(£)|, we write 

-j /»oo i /*oo 

Ff (£) = x / Fq(£-s)£Ff,(£ + is)ds + — Fq{s)CF fl (e + i(£ - s))ds. (23) 

J ^/2 471 Jg/ 2 

By Lemma [2. 121 we have 


|£F^(e + is)| < 2\/2 





ip'(/j 2 )(ld 2 — |e + is | 2 ) 
^(/i 2 ) — ^(|e + is| 2 ) 


1 + V>(s 2 ) 


(for the second inequality observe that the expression under the square root is bounded 
by a constant when s < 2/t and by ijj'(fl 2 )(l + s 2 )/(-0(s 2 ) — -0(/i 2 )) when s > 2 p). The 
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right-hand side decreases with s > 0. Hence, 


'«/2 


Tq(£ — s)CFjx{e + is)ds 


< 


C 3 (e,jd 


(1 + e/ 2 )(l +^(e 2 /4))V 2 


I - s)\ds 


< 


C 3 (e, fi,'i/j)C 2 (£,a,b) 


ds < 


8 C' 3 (e, jl, ip)C 2 (£, a, b) 


(i+$/2)(i+ ^(e/w /2 h/ 2 a + le - *i) 3 - a+o(i+^ 2 )) l/2 ’ 


in the last inequality we used the fact that 4^/>(£ 2 /4) > ^/>(£ 2 ) and that the integral is 
bounded by 1. The estimate of the other integral in (1231) is simpler: | CF^e + zs)| < 
C±{s,p) for all s G R, and hence 


/ Fq{s)CF^{e - s))ds 
U/2 

Therefore, for £ > 0, 

1-77(01 < C 5 (e,a,b,ji) 


^ rit ~\ ^C A {e,iL)C 2 {e,a,b) 

< C i( e,») |^(«)|<fa < 2(1+;/2)2 


+ 


(i + i£i)(i+^(a ) 1/2 (i + iei ) 2 


Since Ff(— £) = the above estimate extends to all £ G R. We conclude that for 

all £ G R, 


(l + ^(£ 2 ))|^/(£)r< 2 (C 5 (£,a, 6 ,/i )) 2 


i + i+i/;(e) 


.(i + ICI) 2 (i + IO )V’ 

and the right-hand side is integrable because (1 + |£|) _2 (1 + - 0 (£ 2 )) is bounded. 


□ 


3.3. Estimates for approximate eigenfunctions. Following [16], we introduce the 
following notation: 

px 

vo(x) = c + / z 2 v(z)dz, 

Jo 

roo 

In= G^(x)dx, 

Jo 

We recall two fundamental estimates, which were proved in [TB] for -0(£) = (£ + l ) 1 / 2 — 
1, but their proofs work for general non-constant complete Bernstein functions if} such 
that V’(O) = 0. One minor change is required in the proof of Lemma 13.41 an extra 
term M 2 c appears when Lemma 13.21 is applied (as compared to the application of [ I6l 
Proposition 4.2] in the proof of | H6l Lemma 4.2]). This extra term is absorbed into 
M 2 u 0 (b). Also, note two typos in the first displayed formula in the original statement 
of |161 Lemma 4.2]: the norm in the left-hand side should not be squared, and the term 
\ n Gji n: b(a) is missing in the right-hand side. (These typos did not appear in the the other 
displayed formula in the original statement, which was the one used later in the proof of 
the main result.) 

Lemma 3.4 ( jl6l Lemma 4.2]). We have 

W^-dFu — ^n<Pn\\L 2 (D) < C(a, b, -0)(( 1 + A n)Gp, nt b(a) — G'^ nb {a) + G'^ n b (a) + Ip, n + A-). 


Voo(x) = / v{z)dz , 


G>,b(^) — Gfj,(x — b) + G^(x + b). 
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More precisely, we have 

II^D^n II L 2 (D) — 2(o b) 


G^(aW 26) + p(2b)I ^ + 2-(a) 


26 2 


l^r 


2b 


(G^ b (a) - 2bGh n b {a) + b 2 G"^ b {a))v,{b) 
b 2 

I O/O t \ fu\ I /r I p,n,b( a ) , 2z/(a) 

+ 2Gp, n 'b(a)voo(b) + v(b)Ijj, n H---h 


l^n 


Lemma 3.5 ( [151 Lemma 4.3]). VFe have 

|II^IIl 2 (d) — a | — 8(Ip, n 4" 1/ pn)- 


More precisely, 


a — 


Sinfe) - 4/ s „ < fc||| 1(D) < a + ^ + 4f fc (l + sin^J. 


I^n l^n 

Lemma 3.6. If if is unbounded, then for n > 2, 

C(a , b, if) 


II^Dy^n ^nPn | L 2 {D) — 


n 


and 


20a .. Il9 

a-< || (p. 112 


n7i 


n|lL 2 (D) — 


^ a -1~ 


36a 


n7r 


Proof. By [T9], Lemma 4.21], 


COS'#,, 


4 = £G>(0) =-^ - £F„(0 + ) = 


COS I? ^ hf'(p 2 ) . 1 


< - 


1f(p 2 ) p 


p p 

Furthermore, by complete monotonicity, 

4 > f G^{z)dz > f {G^{x) - G'^(x)(x - z)+ \G"^x){x ~ z) 2 )dz 
Jo Jo 

= xG^x) - \x 2 G'^x) + ix 3 G"(x), 


so that 


G^x) < 


px 


G'(x) < 


px 


2 ’ 


< 


//x 


3 ' 


By Lemma 12.131 for p > p 2 , 


if(p 2 )G^x) < 


C(if,x) 


p 


(24) 


Finally, p n > ^ for n > 2. The result follows from Lemmas 13.41 and 13.51 □ 


3.4. Distance to nearest eigenvalue. Let a {Ad) denote the spectrum of Ad- Recall 
that the spectrum of Ad is purely discrete (see Subsection 12.51) . and the eigenvalues of Ad 
are denoted by X n . The following result was given in |16] for if(f) = (£ + 1) 1//2 — 1 only, 
but the proof extends to arbitrary self-adjoint operators Ad that preserve the spaces of 
even and odd functions. 
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Lemma 3.7 ( [[T5I Proposition 4.2]). We have 


dist(A n , cr(A D )) < 


ll^n||L2(D) 


( 25 ) 


In fact, if A e f) en and A°f) d are the restrictions of Ap to the (invariant) subspaces of L 2 (D) 
consisting of even and odd functions, respectively, then (1251) holds with ct(Ad) replaced by 
cr(A e f) en ) when n is odd, and by a(A 0 jf d ) when n is even. 


The following result is an immediate consequence of Lemmas 13.61 and 13.71 


Corollary 3.8. If if is unbounded, for all n > 7 there is a positive integer k{n ) such that 

C(a, b, if) 


| A n A k{n) | ^ 


n 


Lemma 3.9. Suppose thatlim^oo^if'^) = oo. Forn larger than some (integer) constant 
C(a, b, if) the numbers k(n) are distinct. Moreover, for any e > 0, for n larger than some 
(integer) constant C(a, b, if ,e), 

if{(j2 n - e) 2 ) < A fc(n) < if ((fin + e) 2 ). (26) 


Proof. Let £ G (0, ^). For some ( n G (/f n ,/J n + e), 

if((/u n + e) 2 ) - if(/u 2 ) = 2ei n if\e n ). 
Since < ff 1 + £ < —, it follows that 


\2\ „/./ ~2 \ ^ 


if((ji n + Ey) -if(fi n ) > 


im 


Since ^ 2 1 ^ , we have lim^oo f 2 if(f 2 ) = oo, and so, by Corollary 13.81 for n greater 
than some constant C(a,b,if,e), 

yiifpn T £) ) VKhn) ^ l^n Afc( n )|. 

Since if> is concave, 

y{pl) - y((Pn - ef) > lf{{Pn + e) 2 ) - lf(j%). 


Finally, A n = if(jln). This proves (l26lh 
Observe that, by Lemma [2. 101 

7 T 7T 3 . ~ \ tt 6 a . _ . 

a Fn +1 a pn ^ d" np. n np.n+1 — ~2 (/hi+1 ^ n — "2 Tfi Ij ^^+1 ’ 

so that fi n+ 1 — fi n > ^(1 + ( TO _ 6 1 ) 7r )^ 1 > for n > 3. The first statement of the lemma 
follows hence from (12^1 ) with £ = JC □ 

Lemma 3.10. Suppose that lim^,*, (if'if) = 00 . Then k{n) > n for infinitely many n. 


Proof. By Lemma [3.91 

A k(n) P yiihn — 16 ^)") 

for n large enough. On the other hand, by (fT2]h 

A„-i < 

for all n > 1. Finally, by Lemma [2.91 and Lemma [2.101 I)f ln < ^ ^ for infinitely many 

n, and hence 

~ _ 7T _ rm _ 1.Q _ 7r ^ mr _ /37 t _j_ n \ _ n _ (n—l)n 

n 16a 2 a a i^ n 16a 2a ^ 8a 16a ' 16a 2a 

for inhnitely many n. □ 
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3.5. Trace estimate. Recall that the kernel functions of the operators exp(— tA) and 
exp(— tA D ) are denoted by T(t]x — y) and T D (t]x,y), respectively. Furthermore, 0 < 
T D (t ; x, y) < T(t;x — y) for all t > 0 and x,y E D = (—a, a), and the Fourier transform of 
T(t] x) is exp(—+(+)). In order to estimate the number of eigenvalues A n not counted as 
A k(n) for n large enough, we use the trace estimate method, applied previously in in [18[ 
Section 9], [201 Section 5] and [IB; Step 4 of the proof of Lemma 4.4], see also [ 2 j [17 ] . 

Lemma 3.11. Suppose that lim^oo £+(£) = oo. For n greater than some constant 
C{a, b, 0) we have k{n ) = n. 

Proof. Let e = ^ and let N be the constant C(a, b, in Lemma 13.91 Dehne J = 
{ k(n ) : n > N} and let J’ = {j > 1 : j J}. We claim that it suffices to show that 
| J'| < N. Indeed, there is no > N such that k(n o) = 1 + max J', and k(n ) is strictly 
increasing for n > N. It follows that k(n) = k(n 0 ) + n — n 0 for n > Uq. If | J'\ < N, then 
kfrio) = \J'\ + (no — N) < n 0 , so that k{n) < n for n > n 0 . Since k{n) > n inffiiitely 
many times by Lemma [3. 101 necessarily k{n ) = n for n > no, as desired. 

Let t > 0. By the assumption, 0(£) > jl°g£ — C(t) for some constant C(t), and there¬ 
fore exp(—b0(£ 2 )) is integrable. Therefore, T(t\x) is bounded in x E R. In particular, 
Td(£; x, •) is in L 2 {D), and so, by Parseval’s identity, 



(T D (t; x, y)) 2 dydx 



T D (t', x, y)ipj(y)dy 


dx 



e 2Xjt (ipj(x)) 2 dx = 2Xjt . 

3 = 1 


On the other hand, by Plancherel’s identity, 


/ a pa poo r\ poo 

/ {T D (t-,x,y)) 2 dydx < 2a / ( T(t; x - y)) 2 dy =— e~ 2t ^ 2) df. 

-a J—a J—oo ^ JO 


It follows that for all t > 0, 


oo 


£ 


\ t 2 a 

e~ Xjt < — 

7T 



(27) 


Observe that 


^ ] g AP — ^ ] g ^k(n)t 


> g-b((/In+d 2 )i > 


) -i/>((n7r/(2a)+e) 2 )t 


j'E J n=N n=N-\-l n=N 

Denote = mi/{2a ) + £ = (n + 4 ) 7 r/( 2 a). Since eis concave in z > 0, 


rO+i 




fO+i 


y2 _ y2 
Sn+1 S ■ 

72 ITT 

Sn+1 Sn 


e + J-% e -^(£ 2 + i) ) d£ 


2^+i-Un+l-^n - 


3(£n + £ra+l) 


II e -hH£ 2 ) _|_ 


F 2 - F 

Sn+l S n 

£n+l £ ^n^n+l 2,f n _ 


3 (£n + £n+l) 


n P -t'KSl+i) 


= — ( 3n + 3 e -mi) _ 3n + 2 e -^« 2 +i) 
2a V 6 n + 5 6 n + 5 
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2 a 


7 r 




'bv 


V < E 


< 


n=N 

3N + 3 


+ 6 n + 5 6 


6AT + 5 

3iV + 3 


3n + 3 e -^(d) _u 3n + 2 p -m 2 n + l) 

6n + 5 

oo 

e -m 2 N ) + 


i-wa 2 ) 


< ° iV e -mi) , v 
“ 6 AT + 5 ^ 


n=JV+l 


j&J 


(the second inequality is a consequence of < 1 , while the last one follows 

from A fe(n) < + e) 2 ) < ip(£, 2 ) for n > N ). By Q2Z}, 


J2 e~ Xjt < — f e~ m2) d£ - J2 e ~ Xjt < — e-W&dt + 

jeJ' ^ jeJ n 


3N + 3 
6 iV + 5 


Passing to a limit as t —> 0 + , we obtain 


, ,, 2a 3iV + 3 Ar 1 3AT + 3 Ar i 

\J \ E -6v + , r ~ ^ + o + a at , r < ^ + 1- 

7T DiV + 5 3 DiV + 5 


This shows that | J'\ < N, as desired. 


□ 


Proof of Theorem \l.l[ By Lemma [3.111 k(ri) = n for n large enough. Hence, by Corol¬ 
lary [XS 1 


An — An + O(^) — ^(An) + ^(^). 


□ 


3.6. Properties of eigenfunctions. As in the previous articles [16] [18] 20], the inter¬ 
mediate results in the proof of Theorem 11.11 provide some approximation results for the 
eigenfunctions. The details of the argument differ slightly from that of ns HOT , so we 
sketch the proofs. 

Proposition 3.12 (see [201 Proposition 1 ] and [ IB] Proposition 4.9]). Suppose that 


lim 


£—>•00 


£^'(£) = oo. With the appropriate choice of the signs of <p n and with 


fin || Pn || L 2 (D) 

we have fi n = \fa + 0 (£) as n —>• oo, and 

1 


|| Pn fin Pn ||l 2 (D) O 


irW(ir) 2 ) 


as n —>■ oo. 


Proof. By Lemma [3761 indeed fi n = y/a + 0{~). Let a n j = (<p n , (fij)L 2 (D), so that <p n = 
YlgL i a n,jPj in L\D). We choose the sign of tp n so that a n ,n > 0. We have 

\\Pn finPn | £2 E ||^n 0^n,nPn\\L 2 (D) T (o^n fin \ 

II An &n,nPn || L 2 {D) T 11| 0^n,n,^P n || L 2 (D) || An||L 2 (D) I 

E 2||<^n 0^n,nPn | L 2 {D) ■ 

As in the proof of Lemma [3.91 for n larger than some constant, if j 7 ^ n and s — then 
|Aj — A n | > max('0((/r n+ i - e) 2 ) - + e) 2 ), A((An - e) 2 ) - if ((jin- 1 + e) 2 ) 

> 2 <=iHV((<=£*) 2 ) • (£ - 2e) > £ f <A(f ) 2 ). 
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Therefore, 


II Tn 0i n,nPn\\L 2 (D) ^ ' l^re,. 


< 


c 1 


j¥=n 




2a 


£ a - t>) : 

j¥=n 


a. 


n,] I 


< 

— mr 

2 a 


Cl 


C 2 {a,b,ip) 


-0'((1HE)2) A «^||l 2 (D) — ^nTr^^/^nTT ’ 


again by Lemma 13.6 


□ 


Proposition 3.13 (see [2D' Proposition 1] and [ T51 Proposition 1.2]). Suppose that 
linic^oo £0' (£) = oo. With the appropriate choice of the signs of ip n and with 

(—l)^" 1 )/ 2 -^ cos (jx n x) when n is odd, 
fn{ x ) = { / 2 j_ 




when n is even, 


we have 


II fn ~ || L 2 (D) — O I —j= + 


n 


n-K\2^l((rm\2\ 


as n — y oo. 


Proof. Clearly, 

II fn ~~ Pn\\L 2 (D) C || fn ~ n\L 2 (D) + (yjll^n — PnTn ||l 2 (D) + | — 111| Pn || L 2 (D) ■ 

The middle summand is C(l/((^) 2, 0'((^) 2 ))), while the last one is 0(4). Finally, by 
the definition (1221) of (p n and the properties of q{x) and F lt (x), 

\W~afn ^n|li 2 (D) — / (q(~ x )Gp n (a + x)- (-1 ) n q(x)G firi (a - x)) 2 dx 

J —a 

roo /•oo 

<4/ (G~„ n (s)) 2 ds< 40^(0) / G-, n {s)ds = 4G^)CG fl M. 

Jo Jo 

Since 0^(0) = cost? a < 1 and £G M (0) = < 4 (see (1241) h we have 

WVafn - ^n|U 2 (B) = 0{^=). □ 

Proposition 3.14 (see 1201 Proposition 2] and [16] Proposition 1.3]). Suppose that if 
£2 > Ci > 1, then 


> M (— 

m) ~ Ui 

for some M,e > 0. Suppose in addition that 

lirninf £ 3//4 0'(£) > 0. 

£—>■00 

Then p n {x) are bounded uniformly in n > 1 and x 6 (—a, a) 


(28) 


(29) 


Condition ( 1251 ) is known under various names, including weak lower scaling condition 
and subregularity, such a function 0 is also said to have positive lower Matuszewska 
index. We remark that although ( 1251 ) does not imply ( 1251 ) . examples of complete Bernstein 
functions which satisfy ( 1251 ) . but not ( 1251 ) . are rather artificial. 

Proof. Observe that £0'(£) diverges to 00 as £ —> 00 , and therefore main results of the 
present article apply. Furthermore, by ( 1251 ) . we have T(t, 0) < Ci(0)-y/0 1 (l/t) for t < 1, 
see, for example, j6j Theorem 21]. 
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We have 

II (D) = e Xnt \\T D (t)ip n \\ L oo^ 

< e Xnt \\ T D (t)(cp n - ^^„)||l°°(D) + e Xnt j^\\T D (t)ip n \\ L o°( D ). 

Since |<^ n (x)| < 2, the latter term in the right-hand side does not exceed j-e Xnt . For the 
former one, observe that \T D (t)f(x)\ < \\T D (t,x,-)\\ L 2 {D) \\f\\ L 2 {D) ,T D (t,x,y) < T(t,x-y ), 
and, by Plancherel’s theorem, 

1 f°° 

iin«,-)iitf ( B, = je-‘* K yd( =rpt, o>. 

Finally, T(2t,0) < C'i(^)y / '0 _1 (l/(2t)) < Ci(ij >)when t < 1. Therefore, with 
t = j~, 

An 

\Wu\\l°°(D) < ^(Ci(l/>)) 1/2 (^~ 1 (A rl ,)) 1/4 ||/3n¥ , n - <Pv\\l 2 (D) + fp 

In the right-hand side, /3 n = 0(1), ^ _1 (A n ) < (^) 2 (by f[T2jl ). and, by Lemma [3761 

II Pn<Pn ~ ^|U 2 (D) = O ‘ ^ 

Appendix 


Proof of Lemma HOI Let x > 0, 0 < e < \C\ and y > 0, and denote for simplicity 
£ = —e~ l£ x. By the representation (J 5 J) of the complete Bernstein function / and Fubini, 
we have 


Im(/(£)c/(-£))/i(x)dx = c / Im(£0(-£))/i(a;)da; + c / Im(p(-£))/i(a:)da: 


^ J (0,oo) J 0 


Im 


&(-£) 

£ + - 


h(x)dx 


m(dz) 


(30) 


(an estimate which allows us to use Fubini is shown below). Our goal is to provide 
estimates for the integrands and find their pointwise limits as e —> 0 + in order to apply 
dominated convergence. 

For the first integral in the right-hand side of f|30jl . we simply use |£g(—£)/r(o;)| — 
xG(x)H(x ), integrability of xG{x)H(x) and Im(£g(—£)) —> 0 as £ —» 0 + . By dominated 
convergence, the limit as e —> 0 + of the first integral in the right-hand side of (13011 is 
zero. Similarly, \g(— £)/i(£)| < G(x)H(x), G(x)H(x) is integrable and Im(g(—£)) —* 0 as 
£ —y CO, and so also the second integral in the right-hand side of (j30j) converges to zero 
as £ —* 0 + . 

To estimate the last integral in the right-hand side of (001) . we consider separately two 
cases. When x < | or x > 2y, we have 


£ 

£ + 1 / 


<?(-£) 


< T —-—7 xG(x) < —-— xG(x ) 

\x — y\ x + y 

< 3min(l, xy~ 1 )G(x) < 3min(l,y _1 )(l + x)G(x), 


so that by dominated convergence, 


f-y/2 />oo N 

L J2y 


Im 


£ + 1 / 


s(-£) H x ) 


dx < 3 C 3 min(l, y ) , 


lim 

£—>■ 0 + 


ry /2 |*oo N 

/ + / 

>0 J2y 


Im 


£ + y 


( 31 ) 


y(-£) h(x)dx = 0. 
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When f < x < 2 y, we need a more careful estimate. Observe that 

&(-£) _ yg(y) - (-Qy(-Q _ yg(y) 

£ + y v - (-0 £ - y' 

The estimate for g and Cauchy’s integral formula for g' easily give 

\g'(z)\ < C^y~ 1 G(y) 

in {z G C : | Argz| < \C\, y/2 < \z\ < 2 y}, with C 4 = 4Cf 1 . By the mean value 
theorem, 

yy(y) - 


when | < x < 2y, and therefore, by dominated convergence, 

'yg(y) - (~0s(-0 


2 

,2y 


0/2 


Im 


y - (-0 


h(x) 


dx < |C' 4 y 1 G(y)H(y) < |C' 2 C' 4 min(l ,y x ), 


r^y 

lim / Im 


£->0 + 


0/2 


( = 0 . 
V 2/ — (—0 / 


(32) 


Finally, if P t (s) and Qt(s) denote the (classical) Poisson and conjugate Poisson kernels 
for the half-plane, then 

1 


Im 


= 7T cos(e)P ysin£ (x -y cose) + tt sm(£)Q ysilie {x -y cose). 


Clearly, P ys in£ (x — y cose) 1(^/2, 2 y)(x)dx converges weakly to 5 y (x), and therefore 
r?y 

\tt cos(e)P ysin£ (x - y cos e)yg(y)h(x)\ dx < tt yg(y)H(y) < C 2 tt min(l, i/ _1 ), 


0/2 


-2 y 


lim / 7rcos(e)P ysine (x - y cose)yg(y)h(x)dx = nyg(y)h(y). 


£-S>0+ 


0/2 


Furthermore, |tQ*(s)| < ^ and tQ t (s ) —> 0 as t —* 0 + , and hence, by dominated conver¬ 
gence, 

r^y 

\tt sm(£)Q ysin£ (x - y cose)yg(y)h(x)\dx < § yg{y)H(y) < \C 2 min(l, jp 1 ), 


0/2 


£->■ 0 + 

We have thus proved that 


r2 y 

lim / 7rsin(£)Qj /sine (x — y cos e)yg{y)h{x)dx — 0. 

S //2 


/•2y 

0/2 


in. 1 -Mb 1 /.(x; 


£ + y 

r^y 

lim / Im 

e_5-0+ J1//2 


dx < < 72 ( 71 + f)min( 1,t/ ), 


yy(y) 

£ + y 


(33) 


h(x)dx = iryg(y)h(y). 


Due to estimates (13T]) . (132|) and (T33|) . as well as the integrability condition on m, indeed 
we could use Fubini in (EDI)- The same estimates allow us to use dominated convergence 
in the limit as £ —> 0 + . We conclude that 


lim / Im(/(£)y(-0)Mx)dx = tt 

s^0+ ./n 


( 0 ,oo) Jo 


g(y)h(y)m(dy). 
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This proves the first equality in (RI]) . The other one follows by replacing the pair g(z), 
h{x) with 1 and g(x)h(x). □ 
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